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Abstract. The COBE and the first-year WMAP data both find the CMB quadrupole 
and octopole to be anomalously low. Here it is shown, that a finite, multi-connected 
universe may explain this anomaly, supporting earlier analyses 0 m A novel 
technique, pioneered by m is used to compute the spectrum and its variance up 
to k=102. Based on the properties of the Lie group of rotations of it is shown that 
the spectrum and its variance may be computed solely from the matrix elements of the 
group-averaging operator, for each of the manifolds S^/I*, /O* and S^/T*. Further, 

it is proved that the spectrum may be calculated solely from the radial function, due 
to the symmetry properties of the Lie-algebra, which is rigorously proven. It is shown, 
that if the topology of the universe is S^/I* the uncertainty on the estimates for ^ltot 
may be improved by an order of magnitude. Finally, the paper highlights how the 
unavailability of an explicit probability function for the observations, given the model, 
is a challenge for Monte-Carlo simulations of the spaces A^/E which has to be addressed 
in future work. 


E-mail: jguSmst. dk 

1. Introduction 

Both the COBE data and the hrst-year WMAP data hnd the CMB quadrupole and 
octopole to be anomalously low mq. A hnite, multi-connected universe may explain 
this anomaly; see jl] for a review. While studies of flat spaces such as the 3-torus 
have not provided a good £t to the observed CMB power spectrum, a preliminary 
study of the so-called binary polyhedral spaces^ showed an excellent £t, particularly for 
the dodecahedral space [S]. Unfortunately that preliminary study suffered two major 
weaknesses: (1) because of computational limitations it computed only the hrst two 
terms C 2 and C 3 of the power spectrum, and (2) it neglected to compute the variances 
of the predicted Ci, making a proper statistical analysis of the results impossible. The 
present article resolves both those problems. A novel approach to the computation 
yields reliable predictions for the power spectrum from C 2 through Cis, along with their 
respective variances. 

f See ini for an elementary introduction to spherical spaces 5'^/r, including the binary tetrahedral, 
binary octahedral and binary icosahedral spaces. 
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For each binary polyhedral space 5'^/r, the computation hnds the eigenmodes of the 
Laplacian on that are invariant under the action of the group F. These F-invariant 
eigenmodes dehne the basic modes of the primordial density fluctuations, before the 
decoupling of the CMB. Exploiting the symmetry properties of these eigenmodes, 
expressions for the Ci and the variance of the Ci are derived, that can be computed 
solely from the radial functions and the matrix elements ( kirn \ G \ ki'm' ) of the 
group-averaging operator G, which projects the space of eigenfunctions on 5^ down to 
the eigenspace of S'^/F. For the binary polyhedral spaces S^/T*, S^/0* and S^/I*, the 
author has carried out the computations up to A; = 100 in the space of 3-dimensional 
modes, yielding the estimates for the Ci and their variances up to £ = 15. 

The ultimate goal of this work is, of course, to compare predictions to observations. 
Here two issues arise. The hrst is the question of what spectrum to use for the initial 
density fluctuations. The present study adopts the neutral choice of a scale-invariant 
initial spectrum, keeping in mind that the justihcation for such a choice is tenuous in a 
small universe and may be subject to future revision. The second issue is the reliability 
of the observed low-£ CMB. Several authors have found unsettling alignments between 
the low-£ modes and the ecliptic plane laiH], suggesting either some hitherto unknown 
solar system effect on the CMB or perhaps some errors in the collection and analysis 
of the data. In light of these uncertainties - and the signihcant revisions they may 
force onto the low-i power spectrum - it seems prudent to wait for the second-year 
WMAP data as well as a plausible resolution of the strange solar system alignments 
before drawing any conclusions about which topologies best £t observations. 


2. Methods 


The space-time metric for the universe in case of both the and the S'^/F models are 
given by the line-element 

ds^ = c^dt^ — R{tydr^ = a{rjf‘{drf‘ — dr‘^) (1) 

where 


dr'^ = dw"^ + dx^ -|- dy‘^ + dz^ 


( 2 ) 


is the spatial distance on the 3-dimensional unit sphere + y'^ + = 1, and 

R{t) = a{ri) is the cosmic scale factor as a function of time respectively conformal time 


V- 

The evolution of the scale factor is given by the matter dominated Friedmann- 
equation (see ffH) 


n 

Ho 


= Qr 


,x 


-1 


- 1 - — iltot + 1 


( 3 ) 


where H = a'/a, x = a/O q is the ratio of the scale factor to the scale factor today and 
7^0 = l^toi ~ Here a' means the derivative of a with respect to the conformal 

time, and fltot = i^mass + is the value of the total energy-density today, expressed 
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relative to the critical density needed to close the nniverse, which has contribntions from 
both matter (baryons and cold dark matter) and from the cosmological term. 

According to the cnrrent nnderstanding, the observed variations in the CMB 
temperatnres arises as a conseqnence of backgronnd flnctnations in the gravitational 
held, which arose dnring the inhationary phase, conpled to density variations in the 
early nniverse. After the nniverse became transparent to electromagnetic radiation, 
the radiation coming to the observer from the last scattering snrface, will have a 
hngerprint from the gravitational held, as well its valne at the point from where it 
was originally emitted, as from the valne of the (changing) held, dnring its travel to 
the observer. Radiation originating at a low gravitational potential <h, will be more 
redshifted from climbing ont of the lower potential, and hence be seen as cooler, than 
radiation originating from a local area with a high potential. And radiation experiencing 
a larger than average dilation ehect from expansion dnring its travel to the observer will 
likewise be seen as more red-shifted than the average. The hrst ehect is known as the 
”Sachs-Wolfe ehect”, the second as the ’’integrated Sachs Wolfe ehect”. Following P 
we get the following temperatnre hnctnations: 

6T{n) = ^^r]Lss, (ho - VLSs)n\ + 2 f — ^^dr] (4) 

o JrjLSS orj 

This eqnation does not inclnde a third ehect, the Doppler-ehect, which arises from the 
relative movement towards the observer, of the matter in the region of last scattering. 

The gravitational held <F here, has its origin in the following ansatz for a pertnrbed 
space-time metrics: 

ds^ = a‘^{ri)[—{l + 2^)drj^ + (1 ~ 2^)'^ijdx^dx^] ( 5 ) 

which is the simplest possible pertnrbation considering only scalar perturbations, and 
working in the longitudinal gauge. On large angular scales, one may ignore the ehects of 
the anisotropic pressures generated by the metric perturbations, and hnd that T = $. 
Likewise, the inhuence of radiation on the metric between the last scattering surface and 
today may be neglected. Under these circumstances, the evolution of <F is determined 
by 

+ mil + cl)^' - + [ 2 n' + (1 + 3cl){n‘^ - i)]<f = o (6) 

The terms with the sound velocity may be neglected on large angular scales, in which 
case we can solve the equation by the following separation ansatz: 

= (7) 

d,s 

where T|^|,(a;) are the eigenmodes of the Laplacian compatible with the topology R^/T 
of the universe, labelled with the index /3, indicating the eigenvalue —(/5^ — 1) of the 
Laplacian, and an arbitrary auxiliary numbering s of the eigenstates. The effect of the 
above various approximations is, that we can ignore the wavenumber dependency of 
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the perturbations, which reduce the dimensionality of the computational problem a lot. 
It is however known, that a wavenumber dependence of about 5% is found in typical 
models of a nearly flat space-time in the wavenumber range employed in this paper (up 
to (3 = 101) [TH] . With the ansatz ©, we hnd that F satishes 0: 

F” + mF' + {2n' + n^ -i)F = t) (8) 

where the derivatives are with respect to the conformal time p. 

In the case where the primordial fluctuations are assumed to be adiabatic, as 
predicted from certain inflationary theories, the primordial held huctuations are assumed 
to be distributed with Gaussian random complex amplitudes, on each of the eigenmodes 
of the Laplacian which are allowed by the topology of the universe, and with a power 
spectrum taken to be scale invariant (equal power per logarithmic wavenumber interval), 
i.e. Hs = 1, or almost scale invariant, — 1 

=FFPiXf,,. PM = ^ ( 9 ) 

Here ap is the overall scale of the huctuations, whereas are random complex 
Gaussian variables with ensemble averages 

13,3^13',s') = dyyiSss' {^/3,s) = = {X(3,sXyi ^s') = 0 ( 10 ) 

Third order moments have vanishing ensemble averages, while the only nonvanishing 


fourth order moment is 

(^/3i,Si^/32,*2^/93,*3^/34,S 4 ) ^l3l03^S\S3dy2(3A^S2SA ^010A ^ S\SA^0203^ S2S3 (11) 

We further expand {x\f3s) = T|^g(x) on the eigenfunctions to the Laplacian on 

{x\(3lm) = = Rf3i{r])Yi^{9, cf) ( 12 ) 

{x\j3s) ='^{x\j3^m){j3^m\j3s) (13) 

im 

and hnd then from (1) the following expression for the temperature huctuations 

6 T{n) = ^ KyeYeUOA){^irn\f^s)Xy^, (14) 

f3sim 

where 

Kyi = (^Ryiiviss) + 2 /" F'{r])Rye{r]o - r])dri\ P^iP) (15) 

By expanding on complex basis functions and using complex random variables, we have 
arrived at a result which gives a complex value for the temperature huctuations. The 
true helds and temperatures is oh course the real part 

26T{e, 0) = E E mKyiYiUO, I f3s)Xys) (16) 

f3s im 
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The a£'m' coefficients are found from this by expanding the temperature on spherical 
harmonics. In the Appendix is shown that utilizing the expressions (uni) for the ensemble 
averages of the random variables, the ensemble average of Ci may be expressed: 


(o)=i: 

m 


1 

2 i + l 


P 


Kl^trace[G0Py{i)] 

2 2£+1 


(17) 


where Gy and Py{t) are projection operators that projects the whole eigenspace to the 
Laplacian belonging to the eigenvalue fd down to the eigenspace of group-symmetrical 
functions, and to the subspace of a particular value of the angular momentum £, 
respectively: 


= i: l/3s> (/3s| pyi) = \l3im) {/3«m| (18) 

s m 

Using that the operator Gy is actually a sum of right Clifford translations, and hence 
commutes with all left-screw transformations ^2] i as well as using the symmetries of the 
Py{i) operator it is proven in the Appendix, that the expression m can be simplihed 
to 

(Ge) = multiplicity(/?) 

p>e 2 

For we get just 

(C,) = E ^ (20) 

p>e ^ 

The significance of (P|) is, that it shows that we do not have to calculate anything but 
the radial function, and its folding in equation to calculate the spectrum, even for 
spaces with nontrivial topologies. The only way the topology makes the result differ 
from the case of is through the second factor in (IT^ . The multiplicity{(3) is the 
dimension of the space of invariant eigenfunctions, respecting the holonomies T for each 
eigenvalue — (/5^ — 1) of the Laplacian, and is known explicitly |llj . 


For the variance of the Gfs it is shown in the Appendix that it is related to the 
matrix elements of the group symmetrizing projection operator Gy as follows: 


Qeu = {GeGe) - (Q){Q.) = - 


2 2£ + 1 2f + 1 




U' 

m,m' 


( 21 ) 


where the matrix M is derived by ’’symmetrizing” the sum of the matrix elements of 
the group averaging operator, as follows: 


^ {l3lm\Gy\(3em') + (-1)^+^ +“ (/3£ - m\Gy\l3e - m') 

2 ^ KyeKyi^ --- 

P ^ 


( 22 ) 


As shown in the Appendix, the matrix elements of Gy can be found solving a certain 
eigenvalue problem based on the fact that the group symmetrical functions are invariant 
to the holonomies of each manifold S'^/F. The matrix elements of the hermitian operator 
Gy are real, and hence symmetric. As another characteristicum, we note that all the 










CMB spectrum for binary polyhedral spaces 


6 


elements of the covariance matrix are non-negative. For Gy is diagonal, and we 
recover the familiar result: 

= {C,C,) - (23) 

In practical terms, we proceed as follows, which is explained further in the Appendix: 
For each grid point (fitoo ^mass) in a 41 by 21 cell grid spanning the range VLtot = 1-004 
to 1.084 and flmass = 0.2 to 0.4, we solve the Friedmann equation Q from x = 1/1085 
(assuming the last scattering occurred at a redshift oi z = 1085) and the equation (jH]) for 
F and F' numerically, which also supplies piss, the conformal time at the last scattering 
surface. Also, for each (3 = 3.. 101 and ^ = 2... 15 we calculate the radial wavefunction 
in 400 points between 0 and tt, using a symbolic calculator to expand the Legendre 
functions. Interpolating F and F', we calculate the integrand in the second factor of 

for the same 400 values of po — V- Doing then the integral by multiplication with a 
precalculated 400 by 400 matrix and interpolating we get both the hrst and the second 
term in 113, for each cosmic grid point. From the Ky£ values we then calculate the 
spectrum from (unD and the variance from m and (1221). 

From the spectrum and the variance we construct an approximate likelihood func¬ 
tion over the {fltot, ^mass) parameter space. The spectrum and the variance is further 
used to calculate the two-point angular correlation function and its variance, for each 
grid point. 

In the Appendices, the methods and mathematical relations exploited in this work, 
are explained further, along with some additional results. 

3. Results 

The best-£t WMAP models are a natural starting point, as they incorporate a lot of 
constraints from the high £ multipoles of the spectrum. As noted in the report on the hrst 
year results from the WMAP team [2] the position of the hrst acoustic peak constrains 
the universe to be nearly hat. However, for models with a nonzero cosmological constant 
there is a geometric degeneracy along the lines of constant conformal distance to the 
last scattering surface in the £2mass £^tot plane, which allows models with topology 
or S'^/P studied in this paper. Actually [2] gives as the best estimate for flfot the value 
1.02 +/- 0.02 which seems to favour a closed universe, but does not exclude either a 
hat universe nor an open universe. The degeneracy means, that the spectrum can be 
equally well htted by assuming a hat as well as a slightly curved universe. 

It is therefore assumed in this paper that the modelled spectrum of the WMAP- 
team coincides with the spectrum for for low Fs, which is in fact verified. If the 
universe has a non-trivial global topology, the high £ behaviour is asymptotically the 
same if the topology is instead S^/T, provided a scaling of |P| is applied to the spectrum 
of the S^/T models. 
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Indeed, scaling the spectrum for calculated by the methods of this paper, to the 
WMAP modelled spectrum (we choose to use the simple pl-model here), shows a very 
good consistency. Applying next a relative scaling of |r| to the spectrum of the S^/T 
models studied here, we can check whether the alternative topologies can provide a bet¬ 
ter £t to the low i behaviour, while still keeping consistency with the WMAP models 
high i behaviour. As shown in Figure |21 the model hts equally well over the different 
values of Qtot so that the ex post likelihood distribution hardly changes the ex ante prior 
used (a Gaussian centred at Qtot = 1.02 with width equal to the standard deviation of 
Qtot of 0.02 reported by the WMAP team). This indicates that the procedure of scaling 
the overall amplitude of the modelled spectrum to the WMAP model makes sense. 
For S^/I* and S^/0*, the ex post likelihood function becomes much sharper than the 
prior distribution, enabling us potentially to use the low multipoles to constrain the 
cosmological parameters more than achieved by the WMAP team, whereas the poten¬ 
tial conclusions in the case of S^/T* are more mixed. The much sharper distribution 
over fltot especially for S^/I*, means that if the issue of the topology is settled, the 
cosmological parameters might be more strictly constrained. For S^/I* the position 
and width of the two peaks in the likelihood occurs at fltot = 1-028 +/- 0.0023 and 
for the smaller peak, at fltot = 1-017 +/- 0.0015, almost the same value found in ng 
by studying the S-statistic introduced in [Sj and explained in the [Appendix K| of this 
paper. As this statistic is heavily biased to explain what we see in the observations for 
the lowest multipole moments, we favour instead the value 1.028 of the right peak, as 
the most likely estimate. 


The conformal distance to the last scattering surface in the maximum likelihood 
peak of the map for S^/I* is found to be r/ = 0.571 which implies that any matching 
circles in the sky would have a radius of 


9 = acos 


\jl- 2^/175 cot(? 7 ) 


59.6° 


(24) 


We could consider the topology to be a discrete cosmological parameter, and hnd the 
relative ex post probabilities for each one, by applying the Baysian principle, as we 
do for the cosmological parameters. Assuming an apriori distribution of equal ex ante 
probability for each topology, the ex post probability, given the data, is very much in 
favour of the non-trivial topologies. In fact we find by such a recipe the following ex 
post probabilities: 


P{syr) = 0.45 P{sy0*) = 0.31 P{SyT*) = 0.13 P{S^) = 0.12 (25) 


This reasoning, however has several weaknesses. First of all, the choice of using equal 
apriori probabilities is a highly subjective choice. And secondly, models of inflation 
suggest an almost, if not completely flat Universe. Nevertheless, it is thought-provoking 
that the likelihood of the S^/P topology has a sharp peak very near the values of 
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Figure 1. The ex post likelihood distribution over Omass (vertical axis) and fltot — 1 
(horizontal axis) for S^/I*, S^/T*, S^/O* and using a Gaussian primer 


^mass, ^tot found from the high i data. 

In the case where we instead apply a uniform primer over the {^tot,^mass)- 

window, we can not assnme that the modelled spectra should be scaled to the WMAP 
model. Instead we have to resort to scaling the spectra to the observed spectrnm (which 
we know will produce bad results for as the low-i multipoles of the observations are 
systematically too low, at least if a £t to higher frmultipoles shall be realised as well). 
We are here bothered by the fact, that we have not in this exercise, for compntational 
reasons, modelled the higher i multipoles, which might more reliably make a fit between 
the model and data meaningfnl (as all topologies should approach the same curve 
asymptotically). Scaling nevertheless mechanically to the observed spectrnm, we arrive 
at the results shown in Figure El 

It is obvious, that the models S^/I*, S^/0* and S^/T* tend to have their opti¬ 
mum likelihood along the well-known geometrical degeneracy line of constant distance 
(in conformal time) to the last scattering snrface. Bnt except for S^/I*, the maximnm 
likelihood regions lye at relatively improbable values of the cosmological parameters, i.e. 
at high flmass or fltot, and even ontside the window in parameter space stndied here. 

The result of this procedure is shown in FigurelUfor the case of the S^/I* topology, 
for the case of the so-called pl-model, which assumes a power law fit to the spectrum, 
using the WMAP best estimates of the flmass and fltot for this model. 
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S3/l‘ S3/T* 




S3/0’ S3 




Figure 2. The likelihood distribution as a function of fltot along the most likely value 
of ^mass = 0.26 for S^/I*, S^/O*, and = 0.25 for S^/T*. Dotted line shows the 
prior distribution, for comparison. Note the different scales on the vertical axes. 



0.0C4 O.OU O.CIJ O.Cji 0.CJ4 0.C54 O.C<4 O.C'.: O.CM 0.0C4 0.014 0.024 0.C34 C.C44 0.C54 C.C<4 O.C'4 0.CS4 

ProbO ProbSS 


Figure 3. The ex post likelihood distribution over flmass (vertical axis) and fltot — 1 
(horizontal axis) for S^/I*, S^/T*, S^/O* and S^, using a uniform primer, and using 
best scaling to the observed multipole moments in the range of = 2 to 15 
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Figure 4. The temperature correlation function in the point of maximum likelihood 
for each of the spaces (compare Figure^]). For S^/I* at flmass = 0.26 and fltot = 1.028, 
for S^/T* at flmass = 0.25 and fltot — 1.044, for S^/O* at ^Imass = 0.26 and 
fltot = 1.044, and for at flmass = 0.26 and Qtot = 1.016. The graph shows 
the modelled ensemble average with its one standard deviation band, along with the 
observed correlation-function, both filtered to a max £ of 15. 


As seen, the angular correlation function is constrained almost to within one stan¬ 
dard deviation, except around 90 degrees and near 180 degrees, which is explained solely 
by the very low quadrupole moment of the observations. Overall, however, the S^/I* 
reproduces the observed correlation function much better than the model. 

This is also shown in Figure El showing the hrst 15 multipole moments, where the 
observed values are almost fairly well constrained to the ensemble average of the model, 
-b/- 1 times the square root of the diagonal elements of the covariance matrix, for the 
S^/I* model, but is far outside 1 standard deviation for the model, for the i = 2 and 
£= 11 moments. 
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S3 I* 


SS T* 


S3 0* 





Figure 5. The spectrum for S^/I*, S^/T* and S^/O* for the optimum likelihood 
values of the cosmological parameters (as in Figure 0J. Fat curves show the spectrum 
+/- one standard deviation (i.e. the square root of the diagonal elements of the 
variance-matrix). Dotted line shows the WMAP model +/- the calculated cosmological 
variance, and thin line shows the model +/- its standard deviation in its optimum 
likelihood point of parameter space. Crosses show observed multipole moments 


Any errors or shortcomings remain however the sole responsibility of the author. 


Appendices 

Appendix A. Rotations on 


Rotations on can be described in complete analogy with usual rotations in E^, by 
considering the rotations to take place in the 4-dimensional Euclidian embedding space 
for S^. Working with the variables n = 0,1, 2, 3 x = (x, t/, z, w) a counterclockwise 
inhnitesimal rotation of the coordinate system in the xy-plane an angle df changes the 
coordinates of x into x' 


/I d(p 0 0\ 
-dcf) 1 0 0 

0 0 10 
VO 0 0 1 y 


(A.l) 


A scalar function 'ffx) transforms hereby into 'fi'x') = '?/’(x), as follows: 


fjfx) = (1 -h z (i0Joi)z/’(x) 


(A.2) 


where Jqi , the generator of the inhnitesimal rotation, is 

d d 
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Considering alternatively to be a function of the coordinates a, 6, 0, 
f X ^ f sin{a)sin{6)cos{<j>) \ 

sin{oi)sin{0)sin{(l)) 
sin{o!)cos{6) 

\ w I \ cos(a) 


(A.4) 


one finds that 


r - 


(A.5) 


A hnite rotation of the coordinate system through the angle 0 transforms ip as follows 
ip'{^) = (A.6) 


If instead we think of examining the value of a function at a position x' that arises by 
rotating the vector x the angle 0, ie 


x = Rxy{4>)y^ 


Rxy 


/ cos{<p) —sin{(p) 0 0 \ 

sin{(p) cos{(p) 0 0 

0 0 10 

VO 0 0 1 y 


(A.7) 


we can do it by rotating the coordinate system the same amount, hnding the transformed 
function, and checking its value at the new coordinate equal to the old, i.e. 

ip{Rxy{ 4 >) ■ x) = e*7'^°i0(x) (A.8) 

If we do two consecutive rotations we have to undo them in reverse order: 


iP{R,MRxy{<P) ■ x) = ■ x) = (A.9) 


Appendix B. The Lie algebra of rotations on 

In the following we state some useful properties of the Lie algebra of rotations on S^. 
The generators Jke = —Jek ,i^k = 0,...,3of the rotations are hermitian and satisfy 
the following commutation relations: 

['Jkii Jim] ~ iJmk (®A) 

for any three different indices k, i, m as well as 

[Jkli Jmn\ = 0 (B-2) 

when k,i,m and n are all different. 


For an observer sitting in (x, y, 2 :, w) = (0, 0, 0,1) the three generators Jqi, J 20 and 
J 12 will be the analogue on for usual rotations around the z-, y and x-axis respectively, 
and the sum of their squares will be the analogue of the familiar total angular momentum 
LJ 


r2 _ t2 I t2 I t2 
~ W <^20 ^ '^12 


(B.3) 
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SRq — Ji2 + J30 SRi — J20 + ^31 SR2 — Jqi + J32 (B.4) 

SLq = J12 — J30 SLi = J20 — J31 SL2 = Jqi — J32 (B-5) 

We note, that the left screw generators here are chosen with the opposite sign than if 
thought of as mowing from {x,y,z,w) = (0,0,0,1) in the direction of the z-axis and 


simultaneously making a left-hand rotation of y towards x. This makes the following 
formulas more symmetric between the SR and SL set. 

The screw generators satisfy the following commutation relations: 

[SRo, SRi] = 2iSR2 [SLo, SLi] = 2iSL2 (B.6) 

and analogous, by cyclic permutation. 

Further, any pair of left and right screw generators commute: 

[SLk,SRe]=0 (B.7) 

The following operators commute with all the generators, and are thus Casimir-operators 
of the Lie group: 

+ Jo2^ + Jo3^ + + ^23^ = 

SR^ = SRq^ + SRf^ + SR^^ (B.9) 

SL^ = + SLf^ + (B.IO) 

Further, we can construct lifting operators 

AR2 = SRo + iSRi SR2-AR2 = AR2-SR2 + 2-AR2 (B.ll) 

AL 2 = SLo + iSLi SL 2 ■ AL 2 = AL 2 ■SL2 + 2- AL 2 (B.12) 


These relations (that holds also by cyclic permutation) show that AR 2 acting on an 
eigenstate to SR 2 with eigenvalue sr is also an eigenstate to SR 2 with eigenvalue sr + 2, 
i.e. AR 2 ’’lifts” the eigenvalue of SR 2 2 units, and analogously for AL 2 acting on an 
eigenstate to SL 2 . In a similar way, the adjoint operators are seen to lower the eigenvalue 
2 units: 

ARI = SRo-iSRi SR2-ARI = ARI-SR2-2-ARI (B.13) 

and analogously for AL^. 

The left screw lifting and lowering operator does not change the eigenvalue of the 
right screw operator and vice versa: 

[ALi, Sflal = [All SB 2 I = \AR2, SLJ = [Afi!,, SLJ = 0 


(B.14) 
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We can thus have states that are simultaneously eigenvectors to SL 2 and SR 2 , with 
eigenvalues in a range slmin, slmin + 2 ... slmax determined by the nonvanishing of the 
result of the lifting and lowering operations. The squared norm of a state resulting 
from application of the lifting operator is (we here use the Dirac-notation explained in 
[Appendix Cj ): 

{si sr\AL\AL 2 \sl sr) = {si sr\SLQ + SL^ + i[S'Lo, ALijls/ sr) = (B.15) 

{si sr\SL^ - - 2SL2 \sI sr) = SL^ - sf - 2sl (B.16) 

where we have used that S'L^ commutes with all the generators and thus is a common 
eigenvalue for all the states. 

Similarly, we hnd that the squared norm after lowering is 

{si sr\AL 2 AL 2 \sl sr) = SLf — sl^ + 2sl (B.17) 

We thus see that the number jdi of different eigenvalues si = —{Pi — l),—{Pi — 1) + 
2,Pi —1 is related to as follows SL"^ = Pp — 1. We would get the same result by 
examining the possible range of eigenvalues for any of the left screw operators, so the 
multiplet of states should also be spanned by a basis of eigenstates of the other two left 
screw generators as well, with Pi different eigenvalues. A similar reasoning shows that 
the range of possible eigenvalues for SR 2 is sr = —{Pr — l),—{Pr — ^)+2,,...,Pr — l 
where Pr is the number of different eigenvalues, related to SR^ as follows SR^ = Pf^ — 1. 

The 2 numbers Pi and Pr completely characterise the multiplet, and we can easily 
construct a matrix representation of the operators, by choosing a basis where SL 2 and 
SR 2 are diagonal. The dimension of the representation is equal to Pi ■ Pr. In such a 
representation the nonvanishing matrix elements can be chosen as 

{si sr\SL 2 \sl sr) = si {si sr\SR 2 \sl sr) = sr (B.18) 

{si + 2 sr\AL 2 \sl sr) = i\jPp — sP — 2sl = al{sl) (B.19) 

{si sr + 2 \AR2 \sI sr) = i\JPr"^ — sr'^ — 2sr = ar{sr) (B.20) 

The above analysis illustrates that any irreducible representation of the commutator- 
algebra of the generators Jm of rotations of is in fact the direct product of represen¬ 
tations of the algebras of generators for rotating E^\ S'0(4) = SO {3) x SO{3), which 
is further exploited in [Appendix Ej (see also m which uses this fact to derive explicit 
analytical solutions for the lowest eigenfunctions of the Poincare dodecahedral space). 
Although we can have abstract multiplets for any Pr,Ph the multiplets corresponding 
to single component wave functions have Pr = Pi = p. One may in fact show, that 
the Laplacian for such wave functions is simply — J^. There are thus P'^ states in the 
multiplet, corresponding to the multiplicity P'^ for eigenstates to the Laplacian on S^, 
with eigenvalue — (/3^ — !)• The usefulness of studying the possible choices of commuting 
operators is that we can use the eigenvalues of the operators to label the degenerate 
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eigenstates of the Laplacian on S^. Where the conventional choice is to use the set of 
operators Jqi, we have found it useful to employ instead the set J^, SR 2 and SL 2 . 


The above explicit expressions ()B.18jl to ()B.20fl determine all the generators Jki- 
There is a choice of relative phases between states involved, in the choice of the lifting 
operators. This is purely conventional, and the expressions may be multiplied by any 
complex number with modulus 1. Our use of the factors of i are convenient, however, 
as then the operators Jqi and J20 which are generators of rotations around the z-axis 
and y-axis respectively are both real matrices: 

SRo = ^{AR 2 + ARl) SLo = ^{AL 2 + All) (B.21) 

SR^ = ^{AR 2 - AR\) SL, = ^XAL 2 - All) (B.22) 

Ji 2 = l{SRo + SLo) J 2 o = \{SRi + SL,) = ^{SR 2 + SL2) {B. 23 ) 
J30 = ~ SLq) J31 = -{SRi — SLi) J32 = -{SR2 — S'L 2 ).(B. 24 ) 

Appendix C. The eigenstates of the Laplacian on 


We can in fact use the rotation operators to calculate the eigenfunctions. We use the 
Dirac notation na, well known from its use in quantum mechanics: 

{f\g) = j f{a,e,4>)g{a, 6, cj))dn (C.l) 

where dfl = sm{6)dad6d(j). This expression is what physicists would call the 

expression for {f\g) in the coordinate representation. If however we have an expansion 
of / and on a complete set of basis functions, call them \(3 i m), we can also calculate 
{f\g) in the following matrix representation: 

I, m) = {(3i m\f) G(A £, m) = {/3 i m\g) (C.2) 


as 


or for short 


(/|9) = i: F{fl,i,,n)G(j3,e,m) 

pirn 


(C.3) 


{ / I 9 > = (C.4) 

Considering | a 0 0 ) to be the name of the delta-function on S^, we can write a scalar 
function ip as 

ip{a, 6,(j)) = {a6(j)\'ifj) (C.5) 

An eigenfunction to the Laplacian, say the 0, £, m state, may then be calculated using 
the rotation operators, as follows 

0,4>) = {a6(j)\(3im) = {O 0 0| e*“'^32g*6'J20g*<^Joi \ f] £ m )(C.6) 
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This expression may most easily be understood, by noting that we simply make rotations 
of the coordinate system in the order 0, 9, a that has the effect that the point {a, 9, 0) 
gets new coordinates (0, 0, 0) 

^ 0, 0) = f){a, 9, 4>) (C.7) 

which is maybe easier to grasp than the alternative 

( 0 0 0 I ^ ^ 0 0 I (-C.s) 

which in conjugated form is 

I a9(f) = Q Q Q ^ 

and which can be seen as the \ a 9 (f) state arising from a rotation of the state | 0 0 0 ) 
in the sequence a, 6*, 0, where there is a minus as we rotate the state, rather than the 
coordinate system. 

How do we calculate the expression ()(k6f) ? Well, the rule, which is so familiar in 
quantum mechanics, is that the only way to calculate the result of a linear operator, 
dehned as a function of another linear operator, is to do it in a basis of eigenfunctions. 
So each of the above rotations must be calculated in a basis of eigenstates for the J 32 , J 20 
and Jqi operators, respectively, making it necessary to transform between these bases. 


These rotations mix the various eigenstates belonging to a beta-subspace, but do 
not mix states with different betas. In the calculation of the eigenfunctions we can 
therefore freely insert projection operators such as or 

for that subspace. In that way we get from equation (I(k6fl (and now suppressing the 
/3-index everywhere) 

9,(j)) = ^ (000|/3 i" m" )( i" m" \slsr){slsr\e’'°‘'^^^\sl'sr') 

,m" ,sl,sr,sr' ,sl' ,i\m' 

x{sl'sr'\i'm'){e'm\e^^^^°e^^^°^ \im) (C.IO) 

Simplifying, we get 

i’yemioi, 9, (j)) = ^(000|/300)(£" = 0 m" = 0|s/ = —sr sr)e^°‘^~ 

sr 

x{sl = —sr sr\i m' = 0 )(£ m' = \ Im) (C.ll) 


We have used the diagonal character of the rotation operators, acting in their respective 
bases, to set indices equal, where appropriate. We have also used, that the only /3, £, m 
combination that gives a '0/3€m(O, 0, 0) 7 ^ 0 is P,i,m = /3, 0, 0. Further, the value is 


V’/ 3 oo( 0 , 0 , 0 ) = (0001/300) = ^ 


27r^ 


(C.12) 


We do not calculate ()(k6f) in the coordinate representation, but instead set up a matrix 
representation, where each state is a vector, and the operators are matrices. 
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As a final twist, we can write as 

yumis^i ^1 0 ) 

= ^ N 0 l{^ = 0 m = 0 \sl = —sr sr)e^°'^'^{si = —sr sr\l m' = 0)Ygm(6', 0) 

sr 

= Ryt{a)Yt^{e,ct)) (C.13) 


where 


— i 


202 


\ 7r(2£+ 1) 


(C.14) 


It is seen, that we have derived a Fourier expansion of the radial function Ryi{a). We 
here used the result for analogous to (I(h6fl . that 

Wm(^,0) = = (OO|£O)(£O|e*®^2°e*'^^“0£m) (C.15) 


where 


(C.16) 


Appendix D. The group symmetrical states 

For the manifolds S'^/F we are studying, there exist a set of Clifford translations, that 
connects any point to each of its |F| ghost-images. These Clifford translations are de- 
hned as the screw-transformations that bring the origin to the images of the origin. The 
Universe may employ either left or right screw translations, which however is immate¬ 
rial for the analysis of this paper. Here and in the following, we settle for right-screw 
transformations. It is convenient to imagine the observer to be situated in the ’’origin” 
{x, y, z, w) = (0, 0, 0,1), which has angular coordinates (a, 9, 0) = (0, 0, 0). 


The coordinate transformation involved for the A; = 0,1, 2, 3 is thus, for a screw- 
translation g which brings the origin to a point with angular coordinates {a, 6, 0) 

X 9 ( 7 ^') Rxy{4^')Rzx{9'jRscrew{(^')Rzx{ 9'jRxy{ 0) X (D.l) 


Here, the rotation matrices R^y and Rzx are the usual 4x4 matrices for rotating a 
vector an angle 0 and 6 in the xy plane and the zx plane respectively, leaving the other 
two coordinates unchanged, while Rscrew performs a right-screw transformation along 
the z-axis: 


R, 


a 


/ cos{a) —sin{a) 
sin{a) cos{a) 
0 0 

VO 0 


0 0 \ 

0 0 

cos{a) sin[a) 
—sin{a) cos{a) J 


(D.2) 


The logic is, that hrst we apply a rotation of the coordinate system (hrst 0, then 6 ) that 
brings our target point to lie on the z-axis, then we do the screw, and transforms back 
to the proper 9, 0 direction. 
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Instead of transforming the coordinates, we use our insights from Appendix A[ 
equations (IA.8f) and (IA.9f) . that we can alternatively hnd the transform of the function 


We can make the average over all the |r| ghost images of x, and hnd that it is 

■pL X!^(5 '(x)) = (D.4) 

dig din 

where the sum is over all the coordinates a, 6*, 0 of the ghost images of the origin. We 
see from this, that any group-symmetrical function fj, for which 'ip{g{'x)) = '0(x) for all 
g, evidently is an eigenfunction to the group-averaging operator, with eigenvalue 1: 

(x|V’) = (x|G/3 |^/>) (D.5) 


where 


- — Joi „ — iOnJ20„iOinSR2„iOnJ20a^4>nJoi 


Gs = t^ 


n 


(D.6) 


It may be shown, that the group averaging operator is in fact a projection operator, i.e. 
an operator with eigenvalues either 0 or 1 m Finding the group symmetrical functions 
thus boils down to hnding the eigenvectors to the matrix Gy, that have eigenvalue 1. 
Denoting these |/9 s) where s = 1,... ,multiplicity{f3), the permissible eigenfunctions 
for the Laplacian on S^/T are 


^p,^x) = {x|/3 s) 


(D.7) 


The fact that Gy is a projection operator means that we can write 

G0 = Y.\P‘>){P‘>\ (D.8) 

S 

The operator G has many symmetries. First of all, it is evident from equation jnn), that 
G is a sum of right-screw Glifford-translations. As any right-screw Glifford translation 
commutes with any left-screw Glifford translation 1121 this has the consequence, that G 
commutes with SL 2 as well as the lifting-operator AL 2 . This implies that eigenstates 
to G may be chosen as simultaneous eigenstates to SL 2 . Each such eigenstate will then 
be a superposition of states with identical left-screw eigenvalue but different right-screw 
eigenvalues: 

|/3s) = \l3sls') = ^aj|/3srs/) (D.9) 

sr 

Here we can choose aj to be independent of si as is easily seen by acting with the 
lifting operator AL 2 on the sum. This is seen to be consistent with the known fact 
that the multiplicity, i.e. the number of eigenstates for each /3, is a multiple of (3. Now 
the group-symmetrical functions must be eigenstates with eigenvalue 1 to each of the 
Glifford-translations in the group F dehning the space /T. Specihcally, if we choose our 
coordinate system to have its z-axis aligned along one of the basic Glifford-translations, 
they must be eigenfunctions with eigenvalue 1 to the screw-translations along the z-axis 
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over an angle ar = ;^ = ^ for the case of F = I*, 0*,T* respectively. 

As each of the sr, si eigenstates is an eigenstate to such translations with eigenvalue 
e* “r we realise that the only admissible sr values in the expansion (ID.9jl are sr = 0 
mod Nr- 

Appendix E. The im states 

We can hnd the simultaneous eigenvectors \j3 i m) of the operators and Jqi (the 
conventional i, m set) by simply calculating the eigenvectors of the matrix + Jqi 
the eigenvalues of which are all different. Our software (Mathcad) then automatically 
supplies real eigenvectors, as the matrix is real, meaning that the transformation matrix 
from the |sr, si) basis to the \f3 i m) basis becomes real. We can fix the sign of each £, m 
eigenvector, by calculating the wave function (nm) in a single point, and comparing 
with the following analytical expression: 


{a6(l)\(3 im) = i"^Ry£{a)Yim{0, 0 ) 


(E.l) 


The factor of will disappear if we choose to omit the i in (27) and (28), which is the 
convention used in the following. 


Here are the usual spherical harmonic functions, with the symmetry 


Wm = (-l)”*W-n^ 


(E.2) 


whereas the radial function Rpi has the explicit analytical expression jH] : 



(E.3) 


where the normalisation factor is 


e 


\ n=0 

and PL is the Legendre function 




Due to the factor of (—1)"' with half-integer n, the Legendre function is purely imaginary. 
Combined with the factor of of the complex conjugate of our radial function is 



This means that also 



We note that an extra factor ofi ^ is used in this paper compared to [H] to get consistency 


with ()C.13fl . 
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Using numerical eigenvector-determination to find the transformation matrices 
{sr sl\£ m) has its limits, due to the huge size of the matrices by matrices). 
Memory constraints thus limit the feasibility to a max (3 of 43, on a PC with 512 Mb 
ram. The £, m states can however easily be expanded analytically on the sr, si states. To 
see this, note that the commutation relations (14) show that considered as 3-vectors both 
JL = ^ and JR = ^ satisfy the usual commutation relations for angular momentum 
operators 

[Jx, Jy] = iJz (E.8) 

Eigenstates of SR 2 and SL 2 are thus states with dehnite values of the z-component of 
JR and JL as well as of their absolute magnitude squared, JR^ = JL^ = |(| -|- 1) 
where k = (3 — 1. 


As the components of the JR and JL commute (equation (jR.7jl we can, from 
simultaneous eigenstates of JR 2 and JL^, JL 2 construct eigenstates of their sum, 
L = JR + JL = (J 12 , J 2 Q-I Joi)-! i-e. eigenstates of dehnite = £{£ -|- 1) and L 2 = m, by 
the rule for vector addition of angular momentum j2()j : 

|£rn) = y, ^)|srs/) (E.9) 

srsl 

This shows that the expansion coefficients between the eigenstates |sr, s/) of SR 2 and 
SL 2 and the eigenstates \i,m) of and L^ = Jqi are just Clebsch-Gordan coefficients, 
which have been worked out once and for all [20] • They are real when the relative phases 
are chosen by omitting the factor i in (27) and (28) which we will assume in the following. 


With real Clebsch-Gordan coefficients, the inverse expansion of (IE.9fl reads: 

|srs/) =^(£,|, |;m, ^,^)|£m) (E.IO) 

Inserting the expansion coefficients cg^^gi = (£, |, |; ^^^, y, y) into the expression 
we arrive at the following expression for the 

e,4>) = Y. 0) (E.ll) 

sr 

Using the inverse expansion (jE.1 Ojl . we then End from (IETTI) the following expression 
for the screw-eigenstates {a9(j)\l3srsl): 

i’n .r .,(a. eA) = Y. cg“.iNi,icC ..MS (E.12) 

Isr' 

where m = 

Our group symmetrical functions can then be expressed using the equation ()D.9fl . 
This enable us to determine the coefficients aj by requiring the resulting function to 
be not only invariant to a right-screw along the z-axis, which was used in [Appendix D 
to realise that the admissible sr values satisfy sr = 0 mod but also to another 
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specific Clifford-translation, which together with the translation along the z-axis spans 
the gronp F. For the gronp I* we choose a translation of cTj* = | in the direction 
(6*r,0) = (acos(^),0) and for the gronp T* a translation of af* = | in the direction 
(6*r, 0) = (acos( —|), 0). For the gronp O* we choose a'p* = | in the direction 
(^r,^) = (acos(-^), 0). This generalises a techniqne originally derived for S^/I* in 
|lfij . and for all three binary spaces in ini. The invariance means that 

E <- E 0) (E.13) 

sr isr' 

must be identically zero for all a. This requires all the coefficients in the Fourier- 
expansion in a to vanish, i.e. 

Ec;;'< = o (E. 14 ) 

sr 

where 

CZ' = Y.':aZ.,Ng,cC - e”r“-')y,„(()r.O) (E,15) 

£ 

We can hence determine the a J as the s' eigenvectors with eigenvalue zero, to the matrix 
C^C. Here sr only runs over sr = 0 mod N-p. As the result is independent of si we can 
choose si = 0 when doing the computation. 

Using dEHI) and (lETl) we get the group symmetrical functions. 

We do not in this paper actually use the functions as such, we only use their ex¬ 
pansion for calculating the matrix elements of the group averaging operator. 


Appendix F. The matrix elements {kim\G\ki'm') 


From the results in the preceding sections we hnd that the matrix elements of the group 
averaging operator are 


{kim\G\ki'm') 

= {k i m|(^ \k sr sl){k sr s/|)(^ \k si s'){k si s'|)(^ \k sr' sl){k sr' sl\)\k £' m') 

srsl s' sr' 

— X! {k £ m\k sr 31 ) 0 "^^{k sr' sl\k £'m') 

s' sl sr sr' 


— ^y2m-sl sm. 


2m—sl ^2m' — sl 


a 


^92m'-sl sl 


s'sl 


(F.l) 


Here the sum over sl should only be extended to values (if any) that satisfy 2m — sl = 
2m' — sl = 0 mod N-p. 
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The starting point is the eqnation (HI in the methods section, giving the temperatnre 
as a fnnction of angle in the sky. The eqnation as given, will normally, with complex 
random Ganssian variables Xps, resnlt in a complex temperatnre signal, which is dne to 
the fact that we work with complex eigenfnnctions of the Laplacian. As the observed 
temperatnre is off course real, only the real part of the expression should be retained. 
Further, we have to take into account the chosen relative phases of our wavefunctions 
{a9(f)\pim), given in equation (IE.211 and (lE.Tfl of [Appendix E| Hence we should write 

26T{e,(f) = Ky,Y,,n{(p,e){(3im\(3s)Xy, + c.c. (G.l) 

pirns 

where c.c. is the complex conjugate of the hrst expression, i.e. by (IE.2fl and (IE.711 

c.c= -KfieYe.m{4>,9){-lY-^{Ps\(3im)Xi, (G.2) 
Expanding on spherical harmonics 

ST{9,(P) = J2a,,nYU<P,0) (G.3) 

im 

We hnd 

2a,m = E - {-lY-^{l3s\f3i - m)^] (G.4) 

hs 

From this we get: 

m 

= E 

l3lSll32S2 

X [{(3iim\(3iSi)X/3^si - {-lY~"^{(diSi\(3ii - m)X/3^si] 

X [(/32S2|/32M^ - i-lY-^W - m\/32S2)Xp,,,] (G.5) 


Taking ensemble averages and summing over m, and using (7) in the methods section 
gives 


= '^2\KyeY{l3im\Ps){Ps\pim) 

m m^s 

= Y.2\Ky,Y{l3lm\G0lm) (G.6) 

mfj 

where Gp is the group averaging operator. Hence we have found the following expression 
for the cosmic ensemble average of the Ci. 


(G ^ 1 ^^^ {^j3iTrt\Gp\j3tTn}j 

\^i/ 2-^ cyn 1-1 o 2 -^ I 2 -^ 

m "I” ^ P ^ 


2 £ + 1 


We show in [Appendix Ijthat this may be further simplihed to 


(o) = i E lA>l 


y multiplicity {l3) 


(G.7) 


(G.8) 
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The significance of this relation, which is stated as a conjecture in HHI, but in this paper 
is proven rigorously, is that it demonstrates that the spectrum can be computed solely 
from the radial wave function. The software used (Mathcad) allows to evaluate this sum 
up to /3=101. This limit arises from the fact that Mathcad’s symbolic engine, which is 
used to expand the Legendre functions (lElI), produces exact rational coefficients with 
nominators and denominators that must be smaller than the largest real which is allowed 
in Mathcad. It is a possibility to be explored, that using instead the Fourier expansion 
fl(hl3jl would allow this limit to be raised. 

Appendix H. Cosmic variance of Ci 
We use ()(L4fl and ()E.7fl to get 


/3iSi/32S2/33S3/94S4 

X ms2\(32im)X^^ - - m\P2S2)Xy,s,] 

X {p,ss\Psi' - m')^] 

X ms\Pd'm')XY:Z - - m'lP.spXp^,,] (H.l) 

When we thereafter take the ensemble average, the random variables only contribute in 
pairs (see (8) in the methods section). Doing the pairing between variables in the hrst 
and second factor, and between the third and fourth factor just result in the product 
16(a£mFFn) ^ wMle the other possible pairings produce 


X [-(^-iY'-^’{p,im\Gy,\Pii' - m') - {-lY-^Wm'\Gp,\P^l - m)] 

X _ ni'\Gy,Wm) - {-lY-^{P2i - m\Gy,Wm)] 

X [{Piim\GyYPii'm') + - m’\Gy,\Pii - m)] 

X [{P2i'm'\GyYP2im) + {-lY^^'-"'{P2i - m\G0YP2i' - m')] (H.2) 

One should note, that the matrix elements of G are real for all the manifolds we study 
in this paper, and hence also symmetric (this may be shown strictly), a fact we do not 
exploit here. Summing over m and m' we realise by substitution of —m' for m' that the 
two sums gives identical contributions. Further, the hrst factor is the complex conjugate 
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of the second (albeit real), with /3i substituted for ( 32 . Our dual result becomes then 
that the covariance of is 

11 1 X—^ I ^ |2 


= {C,C,,) - {Cf){C,) = 


E \Km' 


(H.3) 


2 2£ + 1 2f + 1 , 

where the matrix M is derived by ’’symmetrizing” the sum of the matrix elements of 
the group averaging operator, as follows: 

{jSimlGyljSi'm') + (—— m\Gy\j3i' — m') 




ye 


y 


We note that all the elements of the covariance matrix are non-negative. 


(H.4) 


For the matrix (£ m\G 0 \l'm') = Sei'dmm' so that the variance is just 

2 


= 5. 


£ I' 


21+1 


{Ct) 


(H.5) 


Appendix I. The formula for the spectrum 


In the [Appendix F we claimed that the formula for the Gi of equation (lO.Tfl may be 
simplihed to the expression of equation (nm). 

The starting point is equation (1S21), which may further be simplihed by introducing 
the projection operator on the f-eigenspace P{€) = m){i m \ as 




\Km? 


tia.ce[G yP {€)] 


This, as we demonstrate just below, boils down to 

\K 0 i\^ multiplicity(/3) 


(Q) = E 


y>£ 


(5^ 


( 1 . 1 ) 


( 1 . 2 ) 


This equation is very convenient, as it allows us to calculate the spectrum without 
calculating the eigenfunctions which are symmetrical under the holonomy of the man¬ 
ifold considered. The result holds for all manifolds with holonomies that are Clifford- 
translations. 


For 53 we get just 

(o> = E ^ (1-3) 

y>e ^ 

The steps leading from equation jni to o are based on the observation that both 
G, the group averaging operator and P^, the operator that projects on the eigenspace 
belonging to the eigenvalue £{i + 1) of the angular momentum, has certain symmetries. 
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First we note that G is a sum of right-screw transformations, and hence commutes 
with all left-screw transformations, and therefore also with the generator SL 2 of such 
transformations in the wz,xy direction: 

SL2-G-G-312 = 0 (1.4) 

It then follows easily, that the matrix elements of G between eigenstates of SR 2 and 
SL 2 are diagonal in si, because: 

sl{slsr I G I sl'sr') = {slsr \ SL 2 ■ G \ si'sr') 

= {slsr \ G ■ SL 2 \ sl'sr') = si'{slsr \ G \ si'sr') (1.5) 

and therefore: 

{si sr I G I sl'sr') = {si sr \ G \ si sr')6sisi' (1-6) 

Further, the operator AL 2 for lifting si also commutes with G from which we by 
standard arguments can conclude that the diagonal elements are independent of si. 
By construction: 

AL 2 I si sr) = al{sl) \ si + 2 sr) {si sr \ AL^ = al{sl) {si -|- 2 sr | (1.7) 

Noting that AL^ acting on the state | si + 2sr) is proportional to | si sr) and that, by 
complex conjugation 

{si -|- 2 sr I AL 2 I slsr) = al{sl) {si sr \ ALI | sZ -|- 2 sr) = al{sl) (1.8) 

we see that 

AL 2 \ si + 2 sr) = al{sl) \ slsr) (si -|- 2 sr | AL 2 = a/(s/)(sl sr | (1.9) 
From o and ()1.8|1 follows 

al{sl){sl -|- 2 sr I G I sZ -|- 2 sr) = (sZ -|- 2 sr | G ■ AL 2 | sZ -|- 2 sr) 

= (sZ -|- 2 sr I AL 2 ■ G I sZ -|- 2 sr) = aZ(sZ)(sZ sr | G | sZ sr) (I-10) 

demonstrating that the diagonal matrix elements of G are independent of si. 

As the trace of G is equal to the dimension of the space of non-vanishing symmetrical 
states, i.e. the multiplicity, we hnd that, for each of the (3 values of sZ: 

Y-/ ; \n\ 1 \ ^^11 \n\ 1 \ multiplicity{P) 

^(sZ sr I G I sZ sr) = — ^(sZ sr | G | sZ sr) =- (I-ll) 

sr H slsr F 

For the operator Pi, we hnd, hrstly, that the only nonzero matrix elements are the ones 
satisfying the selection rule sZ -|- sr = si' -|- sr' and secondly, although not so easily, that 
the sum 

^(sr sZ I I sr si) (1-12) 

sl 

is independent of sr and hence for each of the P values of sr must be: 

1 2/ -I- 1 

'^{sr sZ I P^ I sr sl) = — ^(sr sZ | P^ | sr sl) = —-— (1-13) 

sl m slsr m 
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Combining the two equations, we realise that 

trace(GPi) = ^ {si sr \ G \ si'sr') {si'sr’ \ Pi\ si sr) 

slsrsl'sr' 

= ^(s/ sr I G I s/ sr){sl sr | | si sr) 

slsr 

because equation (77) requires si = si' and the selection rule for P^ then sr 
which by the independence of the G term on si allow us to write 

trace(GPz) = '^{sl sr \ G \ si sr) ^(s/ sr | P^ | si sr) 

sr si 

multiplicity {(3) 2£ + 1 

which hnishes our proof, apart from proving the above stated properties of Pf. 

To do that, hrst note that SR2 + SL2 is just 2 * Jqi and hence commutes with Pi. 
Therefore 

{si + sr){sl sr I P^ I si' sr') = {si' + sr'){sl sr | P^ | si' sr') (1-16) 

showing that if {si + sr) 7^ {sV + sr') the matrix element vanishes. 

The further remaining detail, of showing independence of J2si{^'^ si \ Pe \ srsl) of 
sr is accomplished by noting that the operator AL2 + AR2 is just 2 * (Jqi + iJ2o) and 
hence commute with Pi. 

We use this to evaluate 

{si sr+2 I Pi {AL2+AR2) I si sr) = {si sr+2 | {AL2+AR2) Pi \ s/+2 sr)(1.17) 

Noting that we have equations for the action of AR2 analogous to the above for AL2, 
we hnd, by letting AL2 + AR2 act to the right on the ket in the hrst expression, and to 
the left on the bra in the second 

al{sl) (s/ sr + 2 I P£ I s/ + 2 sr) + ar{sr) {si sr + 2 | P^ | s/ sr + 2) 

= al{sl — 2) (s/ — 2 sr + 2 I Pf I s/ sr) + ar{sr) {si sr | P^ | si sr) (1.18) 

Next we sum over s/, and note that in doing that, we can replace the s/ — 2 in the hrst 
term on the right hand side of the equation with si (we can do that, because al{sl) is 
zero when si is the maximum value slmax = /3 — 1). The sum of this term then equals 
the sum of the hrst term on the left hand side, and both sums can be neglected. We 
are left with the sums over the second terms 

'^ar{sr){sl sr + 2 | P^ | s/ sr + 2) = '^ar{sr){sl sr | P^ | si sr) (1.19) 

sl sl 

which demonstrates that the sum is independent of sr which completes the proof. 


( 1 . 14 ) 

= sr' , 

( 1 . 15 ) 
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Central to extracting uncertainty-bonnds for the cosmological parameters from the ob¬ 
served properties of the CMB is the setting np of a likelihood fnnction. In the case of 
flat space, or models, this is relatively straightforward, as the a^m’s can be consid¬ 
ered to be independently distributed random Gaussian variables, so that the Cfs have 
a chi-square distribution with 2 i + 1 degrees of freedom. 

In the case of a nontrivial global topology for the universe, as the S^/I*, S^/ 0 * 
and S^/T* studied here, the situation is more complicated. 

The Cfs are then definitely not chi-square distributed, and the Cfs are not inde¬ 
pendent. 

One may in principle calculate the likelihood distributions from the assumed Gaus¬ 
sian distributions of the random amplitudes X by a large number of simulations, for 
each choice of the cosmological parameters. Such grid-based approaches are very de¬ 
manding in terms of computer-time. 

The strategy used by the WMAP team, for nearly flat models, is instead to sample 
the cosmological parameter-space by setting up Markov chain Monte-Garlo simulations 
pg. That approach is critically dependent on the ability to specify the conditional 
probability distribution for the observations, given the model, which is a simple task 
when the G^’s are chi-square-distributed, but not feasible in the case of S^/T. 

In this preliminary analysis it was chosen to use a crude approximation instead, by 
assuming the Cfs to have the Gaussian distribution of cnj. 

From the ensemble-averages of [Appendix G| and [Appendix H| for the Cfs and their 
covariance, such a multidimensional Gaussian distribution is easily constructed from the 
covariance matrix Q of ()H. 3 [ 1 : 

W(Cf\modd) = ^ 

,/det{27^ 

Here the probability depends on the model parameters through the dependence of 
and on the ’’theory”, i.e. the model parameters. In a more refined analysis, one 
has to consider the variance due to measurement uncertainties, but for the low-£ power, 
which we consider here, the cosmological variance dominates. Hence we ignore the com¬ 
plications of measurement uncertainty in the present analysis. 

Often a prior is applied, i.e. an apriori distribution for the model 

parameters, which may be as simple as to constrain the parameters to a certain window 
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in parameter space. In any case, application of Bayes principle result in the following 
likelihood function for the model parameters, given the observed Cfs: 

W{Cf^^\model) * {model) 


C{model\Cr) = 


(J.2) 


Yjmodei W{Cf’^\model) * {model) 

Two different apriori distributions are used in this paper: a uniform primer over the 
cosmic window studied, and a Gaussian primer with means and standard deviations 
corresponding to the WMAP estimate for the pl-model, flmass = 0.257 + / — .025 and 
the general estimate jS] of Qtot = 1-02 + / — 0.02 . 


It’s a valid question to ask, how good the ansatz cm is. It’s clearly unphysical, in 
as much as the Cfs are inherently positive. It’s also known to be slightly biased ra. 
in the case of where the exact distribution is 

W^’{Cf‘\rnodd) ^llx(Cf‘,Cf ,QS) (J.3) 

£ 

i.e. a product of chi-square distributions with mean and 2i+l = 2 jQfff 

degrees of freedom. However, the distribution (HU has the proper mean values, if we 
neglect that the Cfs, are inherently positive ,for all the manifolds 5'^/r, derivable from 
the theoretical relations of the preceding sections: 

(Cf > = cf 

j/^ohs/^ohs\ j/^ohs\ j/^ohs\ 

/ \W /\^i' / —Wi£' 

{R^} = (Y, (Cf - - cf)} = 14 (J.4) 

l,i' 

where the sum over i and i' runs from 2 to 15. 


Doing a direct simulation of (jG.4jl . for an ensemble of 10.000 Universes, which 
takes about 20 minutes calculation time for each choice of cosmological parameters, the 
resulting distributions for the Cfs can be found, and they are reasonably approximated 
by chi-square distributions, with an effective degree of freedom of D = 2 
with some exceptions, however, as seen in Figure [El As an alternative to cm the 
following 2 probability functions have been tested: 

W^iCf‘\model) = n x{Ct\Cf, Q“) (J.5) 

I 

which means that the off-diagonal elements of are simply neglected, as well as the 
distribution: 

W^{Cf^\model) = detP x{Cf\ Gf, 1) = detP x{P ' Cff P ■ Cf, 1)(J.6) 

i i 

where the matrix P is used to get a unit covariance matrix: 

Cf^ = p . Gf = P-Cf P' ■P = (J.7) 

The function (ELI has the proper mean values (j.l.4jl whereas this is not the case for the 
function (HU). 



CMB spectrum for binary polyhedral spaces 


29 


S3/I* S3/T* 




S3/0’ S3 




Figure Jl. The likelihood distribution as a function of fltot along the most likely value 
of rimass = 0.26 for S^/O*, S^/T*, and = 0.25 for S^/I*, with the probability 
function 031) . Compare with Figure [3 of section |3 . Dotted line shows the prior 
distribution, for comparison. Note the different scales on the vertical axes. 


Figure ITT] and Figure QJ] illustrate, that the choice of the probability function W 
has an impact on the calculated ex post likelihood surfaces. Using (l,T.5ll or (l,T.6f) rather 
than (HI has f.x. the effect, that the secondary optimnm for S^/I* near Qtot = 1-017 
becomes somewhat larger, bnt still smaller, than the peak near Qtot = 1-030. Also the 
peak at 1.03 shifts slightly towards lower valnes of Qtot- Hence, the issue of establishing 
a proper probability fnnction W remains an important qnestion for examining the best 
£t estimates for the cosmological parameters in case of the nontrivial topologies S^/T. 

Neither of the proposed distribntions W are exact. Hence the proper solntion in 
fntnre work could be to use that the a^m’s do in fact have a mnltidimensional Gaussian 
distribution, which can be expressed by the covariance matrix Qimi'm' = {demO-i'm')■ 
The inpnt to the likelihood calculation would then need to be the observed a^m’s in¬ 
stead of the observed Cfs. This approach would be complicated by the anisotropy 
of the which are found for all the non-trivial topologies 5'^/F, as the theoretical 
ensemble averages of depend on the alignment of the special directions. If 

these directions can be fonnd by a matching circle’s detection, the procednre conld be 
used quite straightforwardly, however, to get an improved estimate of the cosmological 
parameters. Another intrigning possibility, in case the matched circles search does not 
prodnce a dehnite answer, could be that a systematic rotation in the 3-parameter space 
of rotations on the sky of the a^m’s, althongh compntationally challenging, could yield 
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S3/I* S3/T* 




S3/0* S3 




Figure 32 . The likelihood distribution as a function of fltot along the most likely 
value of ^mass = 0.26 for S^/0*, /T*, and S^/I*, with the probability function 

m- Compare with Figure |21 of section 0. Dotted line shows the prior distribution, 
for comparison. Note the different scales on the vertical axes. 


a best fit alignment between observations and theory, thns identifying any preferred 
cosmic directions, even withont reconrse to finding matching circles. 

It’s rather easy to derive, that lED implies that the last qnantity is distribnted 
as x{R^i 14,14), for all the manifolds S'^/F. This may be compared with the exact dis- 
tribntion of for S^/I* derivable by simnlation from (l(F4ll . shown in Figure [Til 

It is seen that snch a chi-sqnare distribntion does not reprodnce the simulated dis¬ 
tribution of R^ fairly well, for the case of S^/I*. For the more general manifolds such 
as S^/r we do not, however, have an explicit, exact expression for the true multidimen¬ 
sional distribution W{C^^^\model) (but it could be simulated by drawing a large number 
of the random variables Xy^g). 

Nevertheless, as the distribution HD does reproduce the mean values (l,T.4fl if we 
can neglect that the Cfs are positive, we will assume that dm is a workable approxi¬ 
mation, at this preliminary state of analysis. 

A signihcant lesson from Figure iJdl is. that the statistics R^, which is the sum of 
the squared deviation of the observed Cfs from their cosmological ensemble averages, 
measured relative to the expected variance, is a fairly broad distribution. As it is further 
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Figure J3. The simulated histograms for /I* of the moments £{£ + 1) * Ct/2'K for 
selected ^-values (calculated for flmass = 0.26 and VLtot = 1-028). The dashed lines 
show chi-square distributions with the same mean and variance as the histograms. The 
simulated distributions are clearly not chi-square, especially for £ = & and possibly for 
£ = 2 . 


found, that the observed value of B? is fairly close to its expected value of 14, for all 
topologies and across the cosmological parameter space, we can not use as a statistic 
to reject any specihc topology, or any choice of the cosmological parameters. 

Whereas the distribution cnj is very sensitive to deviations of from Cf^ this is 
not the case for the distribution of B^ because of the ’’phase space factor”, i.e. the fact 
that the surface to radius ratio in the 14-dimensional space of the Cfs grows quickly 
with B. This means, that we can use the likelihood distribution (jJ.ljl to pinpoint the 
best values of the cosmological parameters, given the topology, but cannot easily use it 
to discriminate between alternative topologies. 


Appendix K. 2-point angnlar correlation fnnction 

The 2-point angular correlation function C{9) is dehned as the average over the sky of 
C{6) = av{5T{n)5T{n')) with h ■ h' = cos{9), which is related to the moments Ci by 

-1 OO 

C{9) = —Y.{2i+l)CiPeicos{9)) 
e =2 


(K.l) 
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Figure 34. The simulated histogram for /I* of the defined by II.T.411 . calculated 
for flmass = 0-26 and fltot = 1.028 and 10.000 universes. Horizontal axis is R^. The 
simulated distribution is not chi-square, with 14 degrees of freedom, showing that the 
ansatz is not perfect. 


where is the Legendre function. The cosmic expectation value of C{6) is found by 
replacing Ce with {Cf) and similarly, the cosmic variance of C{6) is 

1 CO 

{C{ef)-{C{e)f = j^ ^ {2i + l){2i' + l)Qe,e'Piicos{e))Piicos{e)){K.2) 

As discussed in j3] the observed two-point correlation function is very flat for large angles, 
a feature which the WMAP-models based on flat space or nearly flat space are unable 
to reproduce. The anomaly is mainly a result of the very low observed quadrupole and 
octopole moments. To quantify the anomaly, simulations are reported of the statistic 

fCos{p) 

s{p) = j ^ c{efdcos{e) (K.3) 

It is found, for the best £t AC DM model, for p = 60 degrees that only 0.15% of the 
simulations have lower value than the observed value of S{p). For the running index 
model, the similar result of simulations is found to be 0.3%. 

In this paper, we only study the S-statistic filtered to a maximum value of i of 
15, but the results are quite similar as the S-statistic heavily emphasises the lowest 
multipoles. As an example, it is found that for flmass = 0.26 and Qtot = 1-028 (the 
favoured value from this study) the model with a power law spectrum, gives an 
^(OO) lower than the observed value only in 0.13% of simulations. This is in contrast 
with simulations for S^/P where the simulations performed, for the same values of the 
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Figure Kl. The simulated histogram for S^/I* of the S'(60) statistic, calculated 
for ^Imass = 0.26 and ^Itot = 1-028 and 10.000 universes. Horizontal axis is S'(60) 
times 10“^. The vertical line shows the observed value. The probability of getting 
a value less than the observed is 23%. The simulated distribution is not chi-square, 
but disregarding this one finds the effective degrees of freedom 2mean^/varianse to 
be 2.26, showing that the statistic is sampling very few random variables. 


cosmological parameters, give a 23% chance of a lower S'(60) than observed, see Fignre 

liol 

In jlB] the cosmological expectation valne of the S-statistic is nsed to locate the 
valnes of the cosmological parameters that brings this expectation valne closest to the 
observed valne. 

This might be problematic, however, becanse the choice of statistics is heavily 
biased to emphasize the low valnes of the two hrst mnltipole moments seen in the 
observations, as discnssed in For that reason, it was preferred in this paper, rather 
to nse the maximnm likelihood principle, even thongh it is based on only an approximate 
distribntion W, to search for the most likely valnes in cosmological parameter space for 
each of the topologies. As shown in Fignre ^ and Fignre |21 the maximnm likelihood 
is fonnd for ffmass=0.26 and fltot = 1-028 +/- 0.0023. A secondary bnt smaller peak, 
however, is fonnd near the valne of Qtot fonnd in by minimising the S-statistic. The 
center valne and width of this peak yields fltot = 1-017 +/- .0015. As shown in Figure 
HU and Figure [JU this secondary peak becomes somewhat more pronounced when using 
the alternative approximate distributions dH or (jH). 
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